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Abstract 

This work is devoted to proving existence of global weak solutions for a general 
isothermal model of capillary fluids derived by J.- E Dunn and J. Serrin (1985) [6], 
which can be used as a phase transition model. 

We improve the results of [5 by showing the existence of global weak solution in 
dimension two for initial data in the energy space, close to a stable equilibrium and 
with specific choices on the capillary coefficients. In particular we are interested in 
capillary coefficients approximating a constant capillarity coefficient k. To finish we 
show the existence of global weak solution in dimension one for a specific type of 
capillary coefficients with large initial data in the energy space. 

1 Introduction 

1.1 Derivation of Korteweg model 

We are concerned with compressible fluids endowed with internal capillarity. The model 
we consider originates from the XlXth century work by van der Waals and Korteweg 
[13j and was actually derived in its modern form in the 1980s using the second gradient 
theory, see for instance [121 [2^ . 

Korteweg-type models are based on an extended version of nonequilibrium thermodynam- 
ics, which assumes that the energy of the fluid not only depends on standard variables 
but also on the gradient of the density. Let us now consider a fluid of density p > 0, 
velocity field u E M^, we are now interested in the following compressible capillary fluid 
model, which can be derived from a Cahn-Hilliard like free energy (see the pioneering 
work by J.- E. Dunn and J. Serrin in [6] and also in [HllllH]). The conservation of mass 
and of momentum write: 

d 

—p + dw{pu) = 0, 

Q (1-1) 
— (pu) + div{pu (g)u)- div(2/i(p)D(u)) - V(A(p))divu) + VP(p) = divK, 

where the Korteweg tensor read as following: 

divK = V {pK{p)Ap + 1-{k{p) + pK {p))\V pI"^) - div {K{p)Vp0Vp). (1.2) 
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K is the coefficient of capillarity and is a regular function. The term divK allows to 
describe the variation of density at the interfaces between two phases, generally a mixture 
liquid- vapor. P is a general increasing pressure term. D{u) = ^[Vn +* Vu] being the 
stress tensor, fi and A are the two Lame viscosity coefficients depending on the density 
p) and satisfying: 

/i > and 2fi + NX> 0. 

Here we want to investigate the existence of global weak solution for the system (II. ip . so 
we have to describe precisely the form of the non linear terms coming from the capillary 
tensor. In this goal one can now rewrite K.For this and to simplify the presentation, we 
assume only that k{p) = Kp". With this choice we have to distinguish the case a ^ —2 
and a = —2. We refer to the appendix in section [5] for more details on the formal 
computations on K. We get then: 



i^ij = {AiAB{p) - Al\VA{p)\^ )dij - BadiA{p)djA{p) if a / -2, 

Kij = Ac(Alog(p) + hviog{p)\^)dij - Kdi \og{p)dj log(p) if a = -2. 



with: 



A{p)=p2+\ B{p)=p^+" and Al = -—, AI = / , 



2 + a' " (a + 2)2 ' " {a + 2f 

For the general case, we have A^ {p) = and B' {p) = pK{p) (see the appendix for 
more details). We remark then that the form of the non linear terms appearing in the 
tensor K corresponds to quadratic gradient terms diA[p)djA{p) and some terms in B{p) 
of pressure type. It is the main difference between the system of Korteweg and this one 
of Navier-Stokes compressible where the difficulty state in the treatment of the pressure 
term (see [7j, [15])) Before getting into the heart of mathematical results, one can recall 
first derive the physical energy bounds of the Korteweg system. Let p > be a constant 
reference density, and let 11 be defined by: 

so that P{s) = sn'(s) - n(s) , n'(p) = and: 

dMp) + d\-v{vll{p)) + P{p)d\Y{u) = in p'((0,r) x M^). 

Notice that 11 is convex as far as P is non decreasing (since P' {s) = sn"(s)), which is the 
case for 7-type pressure laws. Multiplying the equation of momentum conservation in 
the system (jl.ip by u and integrating by parts over M^, we obtain the following estimate 
( for more details on the integration by parts, especially on the term divX we refer to 
the appendix section [5]): 



/ [\p\u\' + [Il{p)-Yi{p)) + \^{p)\Vp\^){t)dx+ f f {p{p)\D 

+ e(p)|div^p)dx< / (^ + (n(po)-n(/^)) + ^|V/,oP) 



n)P 



(1.4) 
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with ^(p) = fj.{p) + X{p). We will note in the sequel: 

nt)=[ {lp\u\' + {U{p)-Uip)) + ^\Vp)\')it)dx, (1.5) 
It follows that assuming that the initial total energy is finite: 

£0= [ (^ + (n(po)-n(p-)) + ^|Vpol^)rfx<+oo, 

Jrn 2p 2 ' 

then we have the a priori following bounds: 

n(p)-n(^), and e L1(0,oo,L^(M^)), 

V^^V/9G L°°(0,oo,l2(M^))^, and Vu e L^{Q,oo,M.^f\ 

In the sequel, we aim at solving the problem of global existence of weak solution for 
the system (ILlh so assuming that we dispose from smooth approximates sequel solutions 
{pn,Un)n&i of System (jl.ip . One can remark then easily that the main difficulty lies 
in the fact to be able to pass to the limit in the quadratic term 'SJ A{pn) ® V^(/9„) 
which belongs only to L°°(L^). According to the classical theorems on weak topology, 
V A{pn) V A{pn) converges up to extraction to a measure so how we can prove that 
V = V A{p) VA{p) where p is the limit of the sequence (/>n)neN in appropriate space. 
Notice that if we compare the Korteweg system with compressible Navier-Stokes system, 
the main difficulty for proving global existence of weak solution ( see |15j , ) consists to 
find strong compactness properties on the density p in L^^^ spaces to pass to the limit in 
the pressure term P{p) = ap"' with 7 > > 2 in the case of the works of E. Feireisl. 

In the capillary case, more a priori bounds are available for the density, as V A{p) belongs 
to L°°(i7^(M^)). Hence in our case in dimension N = 2, one can easily pass to the limit 
in the pressure term. 

Another difficulty in compressible fluid mechanics is to deal with the vacuum and we 
will see that this problem does appear in the model of Korteweg, when we want get 
some estimates on Vp. As a matter of fact, the existence of global weak solution for the 
model of Korteweg with constant capillary coefficient for > 2 is still an open problem 
even in the case A^ = 2 if we don't assume some condition on the vacuum as we will 
see in the sequel. The first ones to have studied the problem, are R. Danchin and B. 
Desjardins in [5]. They showed that if we take initial data close to a stable equilibrium 
in the energy space and assume that we control the vacuum and the norm L°° of the 
density p, then we get global weak solution in dimension two. Controlling the vacuum 
amounts here to get a bound of ^ in L°°. Recently D. Bresch, B. Desjardins and C-K. Lin 
in [3] got some global weak solutions for the isotherm Korteweg model with some specific 
viscosity coefficients. Indeed, they assume that p{p) = Cp with C > and A(p) = 0. By 
choosing these specific coefficients they obtain a gain of derivatives on the density p where 
p belongs to L?'{H^). It is easy at that time with this type of estimate on the density p 
to get strong convergence on the term of capillarity. However a new difficulty takes place 
concerning the loss of information on the gradient of u when vacuum appearing (indeed 
the viscosity coefficients are degenerated), so it becomes involved to pass to the limit in 
the term pnUn ® Un- That's why the solutions of D. Bresch, B. Desjardins and C-K. Lin 
require some specific test functions which depend on the density p. This test functions 
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depending of the solution allow to deal with the vacuum. 

Let us mention briefly that the existence of strong solutions for > 2 is known since the 

works by H. Hattori and D. Li [lOllllj . Notice that high order regularity in Sobolev space 

i?^ is required, namely the initial data {po,Pouo) are assumed to belong to x H^~^ 

with s > Y + 4. Moreover they consider convex pressure profiles, which cannot cover the 

case of Van der Waals' equation of state. R. Danchin and B. Desjardins in ^ improve 

this result by working in critical spaces for the scaling of the equations, more precisely 

K K — i E. 

the initial data {po,PoUo) belong to x i?2^i (the fact that is embedded in L°° 

play a crucial role to control the vacuum but or to work in multiplier space). In p[4j, 

M. Kotschote show the existence of strong solution for the isothermal model in bounded 

domain by using Dorea-Venni Theory and TC^ calculus. In [9j, we generalize the results 

of [5] in the case of non isothermal Korteweg system with physical coefficients depending 

on the density and the temperature. We get strong solutions with initial data belonging 

to the critical spaces -62^1 ^ ^21 ^ ^21 when the physical coefficients depend only on 
the density. When they depend on the temperature too, it requires more regular initial 
data to control the temperature. 

Our present result improve the results of R. Danchin and B. Desjardins in [5] and D. 
Bresch, B. Desjardins an C-K. Lin in [3j, indeed we show the existence of global weak 
solution with small initial data in the energy space for specific choices of the capillary 
coefficients and with general viscosity coefficient. Comparing with the results of [3j, we 
get global weak solutions with test function ip G C^(M^) and not depending of the 
density p. Moreover our result can be applied for general viscosity coefficients. Moreover 
in the case of a constant capillary coefficient, we show that we don't need to control 
p in L°° norm to get global weak solution as in In fact we have extracted of the 
structure of capillarity term a new energy inequality using fractionary derivative which 
allows a gain of derivative on the density p. This method shows that the case k{p) = 
K > is a critical case where we can get global weak solutions with small initial data 
in the energy spaces but without any conditions on the density as in [5]. Moreover 
we generalize the existence of global weak strong solutions to the case of capillarity 
coefficients approximating the constant case. To finish we explain why it seems a little 
bit trick to get global weak solution with large initial data in the energy space, indeed 
the problem is then to control some possible concentration effect on |VA(/9)p. 
In section [2] we recall some definitions on the Orlicz space and some classical energy 
inequalities on the system in these spaces. In section [3] we show a theorem of global 
existence of weak solutions in dimension two for specific choices of capillarity coefficient 
with general viscosity coefficients. More precisely we investigate a new structure on 
the capillarity coefficient which allow us to get a gain of derivative on A{p) under the 
condition to get for e > enough small ,,)VA(/?)||j;^cx)(^2) < e uniformly for all x E M 
with r > enough small. By imposing small conditions on the initial data, we obtain 
our results and we explain what happen with large initial data. In the last section we 
investigate the case of the dimension one, and we get a theorem of existence of global 
weak solution in the energy space with large initial data. 
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2 Classical a priori estimates and Orlicz spaces 



2.1 Classical a priori estimates 

We first want recall a priori bounds for initial density in Orlicz space, it means that 
we will work around a constant state p > 0. We rewrite the mass equation by using 
renormalized solutions, and the momentum equation. In the sequel to simplify we let 
P{p) = ap^ with 7 > 1 and a a positive regular function. We get the following formal 
identities: 

JZl^y -P'- ^P'~'^P - P)) + div[n-^(p^ - -p^-^p)] = u . V{p^) 

+ pu • - div(2/i(p)D(u)) • u - V(A(/j))divu) ■n + au-Vp"' ^'^'^^ 

= u ■ divK , 



Notation 1 In the sequel we will note: 

Therefore we find by summing the two equalities of p.6p : 



- div(2p(p)D(u)) • u - V(A(p))divu) • u = u • divK. 

We may then integrate in space the equality ()2.7p and we get (see for more details the 
appendix [5]) : 



/ {p^r + ^-{hip) + \<P)\^P\^)it,x)dx+ldsl 2p{p)\Du\^dx 



X{p)\divu\'^dx < / (po^— H Tj-ripo) + k{po)\V po\'^)ix) dx. 

^R^ JR^ ^ 7 - J- 



(21 



Notation 2 In the sequel we will note: 

|2 



S'^it) = [ (p^ + -^Mp) + lK{p)\Vpf){t,x)dx, 
Jrn^ 2 7-1 2 

^0 = [ {po^-^ + —^hiPo) + 7;K{po)\Vpo\^)dx. 
j^N z 7 — i z 

We now want to estimate this quantity j-y{p) and in this goal we recall some properties 
of Orlicz spaces. 
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2.2 Orlicz spaces 

We begin by describing the OrUcz space in which we will work: 

L^(R^) = {/ G lL(K^)//1{|/|<5} G /1{|/|>5} e 

where (5 is fixed, S > 0. 

First of all, it is not difficult to check that Lp does not depend on the choice of 5 > since 
^ is bounded from above and from below on any interval [(^i, (^2] with < di < 82 < +00. 
In particular we deduce that we have: 

f e l|(M^) if / G L«(M^) and p,q>e. 

e 

Obviously we get meas{|/| >5}< +00 if / G Lp(M^) and thus we have the embedding: 

L5(M^) C L^i(M^) if 1 < ft < g < +00, 1 < p < pi < +00. 

Next, we choose ^ a convex function on [0, +00) which is equal (or equivalent) to for 
X small and to for x large, then we can define the space Lp(M-^) as follows: 

Definition 2.1 We define then the Orlicz space Lp(M^) as follows: 

LliR'') = {fe lL(m^)/*(/) g l1(m^)}. 

We can check that Lp(M^) is a linear vector space. Now we endow Lp(M^) with a norm 
so that Lp{M.^) is a separable Banach space: 

||/||^,(KA^)=inf{t>0/ *({)<!}. 

We recall now some useful properties of Orlicz spaces. 
Proposition 2.1 The following properties hold: 

1. Dual space: If p > I and q > 1 then (L^M^))' = L^,(M^) where q = p = 

p 

p-i- 

2. Ll, = LP + Liifl<q<p< +00. 

3. Composition: Let F be a continuous function on M such that F{0) = 0, F is 
differentiable at and F{t)\t\~^ — ^ a 7^ at t ^ +00. Then if q> 9, 

F(/)g4(M^) ^//GL^(M^). 

Now we can recall a property on the Orlicz space concerning the inequality of energy. 

Proposition 2.2 The function jj{p) is in L^(R^) if and only if p — p e L^. 

Proof: On the set {\p — p\ < 8}, p is bounded from above, since 7 > 1 we thus deduce 
that j'y{p) is equivalent to \p— on the set {\p— p\ < 8}. Next on the set — p| > 8}, 
we observe that for some u G (0, 1) and C G (1, +00), we have: 

uIp-pV <Mp)<c\p-p\\ 

□ 
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Link with our energy estimate 

We recall the definition of the fractional derivative operator A*. 
Definition 2.2 We define the operator A* as follows: A*/ = 

We give now some useful results concerning Sobolev spaces, we start with a proposition 
coming from a theorem of interpolation by Riesz-Thorin. 



Proposition 2.3 The Fourier transform is continuous from in L'^ with p G [1,2] 
q G [2,+oo] and i + i = 1. 

We recall here the definition of homogeneous Sobolev space. 



Definition 2.3 Let s G M. f is in the homogeneous space if: G L^(M^). 

Proposition 2.4 Let f e with s > and f e L^ + L^ with l<p <2. Then f £ L'^. 

Proof: Indeed we have as / G H^: 

[ < +00, 

so /l{|j^>i} G L2(M^). And as / = /i + /a with /i G LP(M^) and G L^. By using 
the Riesz-Thorin theorem, we know that /i G L«(M^) with i + i = l. Asg>2 we then 
have G L^(M^). This concludes the proof. □ 



According to the above theorem and our energy estimate we get that for all T G 



p- p£ L°^{0,T;L^ 



Remark 1 Let k{p) = Kp^ with a ^ —2 We assume now that ^ G L°° if a > —2 and 
p G L°° if a < —2 (These hypothesis will be supposed in the theorem ] 3. We have then 
by using previous properties on Orlicz spaces and 12. 8\) : 

• i/7 > 2 then L](M^) ^ L^(R^) and so p - p £ L'^{H^R^)). 

• ^/ 7 ^ 2 then following the proposition \2.4\ and the fact that L^ = L"^ + L^ we get 
p-pe L'^{H^(R^)). 

We finish with one proposition which give some results of refined Sobolev, for a proof see 

Proposition 2.5 Let 1 < q < p < +oo and a be a positive real number. A constant C 
exists such that: 

II/IIlp < C||/||^-l \\f\f , With P = a{^- 1) and 9 = ^. 
We end with a classical proposition of weak topology. 

Proposition 2.6 Let 1 < r < +oo and K a compact. Let fk a sequel of vector fields in 
such that: 
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• fk is uniformly bounded in L'^{K) and fk converge a.e to f, 
then fk converge weakly to f in L'^(K). 

Proof: We have to show that for all 09 G U' with ^ + ^ = 1 we have: 

r T 

It is enough in fact to show this property for all 99 E ^ with ^ a dense space in U . By 
the Fatou theorem we check that / E L'^{K). Next we set for all A; E N: 

E{k)={yGK: - > 1}, 

and Fp = IJfc>p E{k). We have now = np>o with: 

N = {yGK: yp3k>p \fk{y) - f{y)\ > 1} 

As fk converge a.e to / we have the measure of Lebesgue of N which is null. We choose 
then ^ = {if G U' (K) : suppc/^ c K\N}. We show next easily by a convergence 
dominated that for all 9? E ^, {fk — f, f) -^k^+00 0. And this achieve the proof. □ 
We can now explain what we mean by weak solution of problem (jl.ip in dimension N = 2. 

Definition 2.4 Let the couple {po,uo) satisfy; 

1. po E L](M^), VA{po) E L2(M^) 

2. /)o|noP E Li(M^) 

3. pqUq = whenever x E {po = 0}, 
We have the following definition: 

• A couple {p, u) is called a weak solution of problem on I x with I an 
interval o/M if: 

- p e L°°(L](M^)), VAip) E L°°(L2(M^)), ipp e L\H'^+'' {R^ )) Va e]0,1[, 
and \/ip E C^(IR^). 

- Vn E L'^{L'^{R^)), p\u\^ E L~(Li). 

- Mass equation holds in D' (I x M^). 

- Momentum equation holds in D' (I x R^)^ . 

- limi_o+ /jjjv p{t)f = J^N poV, '^f e ViR^), 

- limi^o+ Irn pu{t) ■ (j) = J^n{pu)o • 0, V(/) E V{R^)^. 

• The quantity is finite and inequality h2. <S'|) holds a.e in I. 
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3 Existence of global weak solutions for N = 2 
3.1 Gain of derivatives in the case N = 2 

We give now a general description of the viscosity coefficients and in the sequel we will 
use this type of viscosity coefficients. 

Definition 3.5 We ask the following properties for the viscosity coefficients A and fi 
which are regular: 

• it exists c > 0, sq > such that Vs such that < s < sq, h{s) > c. 

• it exists ci > 0, m G Z such that Vs > sq, fi{s) < cis™. 

• it exists c' > 0, So > such that Vs such that < s < Sq, A(s) > c. 

• it exists C2 > 0, m G Z such that Vs > Sq, A(s) < C2S™ . 

These hypothesis allow us in fact to control the gradient of the velocity Vn without 
difficulties. Indeed we have then: 

/ / n{p)\Du\'^{t,x)dxdt > c / \Du\'^{t,x)dxdt, 

Jo Jr^ Jo Jm^ 

and we use the fact that: 

/ \Du\^{t,x)dxdt = -{ I / \Vu\^dxdt+ f / \dvvu\^ dxdt) , 
Jrn 2 Jq J^n Jq J^n 

to conclude. In the following theorem we are interested now by getting a gain of derivative 
on the density p. This will enable us to treat in distribution sense the quadratic term 
VA{p)(g)VA{p). 

Theorem 3.1 Let N = 2 and {p,u) be a smooth approximate solution of the system 
with k{p) = Kp°' with a G M and a ^ —2. We assume that if a > —2 then 
i G L°°((0,r) X M^) else a < -2 then p G L°°((0,r) x M^). 

Then there exists a constant ij > depending only on the constant coming from the 
Sobolev embedding such that if: 

I|V/3o||l2(r2) + ||VPo|uo|||l2(R2) + \\h{Po)\\L^ < V 

then we get for all G C§°{R^): 

(j^i+i) < M with < s < 2, 

where M depends only on the initial conditions data, on T, on ip, on s and on ||^||l°° or 
\\p\\l°° ■ 
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Remark 2 When we speak of smooth approximate solution {p,u) of the system ip . 
in fact we can take the solution of 19]. Indeed in this article, we show the existence of 
global strong solution for the nonisothermal system with capillarity coefficient variable in 
function of the density. The initial data are in critical Besov spaces for the scaling of the 
equation. 

In the sequel when I assume "smooth approximate solution (p, u) " I refer to f^. Indeed 
this is compatible with the fact we will need every time of hypothesis of smallness on the 
initial data. 

Remark 3 In fact for the case 7 > —2, we don't need of the hypothesis on the viscosity 
coefficients concerning their behavior around the vacuum as we suppose that ^ G L°° . So 
we can work with degenerate viscosity coefficients. 

Remark 4 We could remark that in the specific case K,{p) = Kp~'^, we get a gain of 
derivative on the density and in particular on log p without condition on the vacuum or 
on the -L°° norm of the density. So the case a = —2 appears as critical. In reality the 
fact that A{p) = B{p) when a = —2 plays a crucial role in the proof. We prove this 
result in the next theorem. 

Remark 5 In fact instead of supposing that ^ G L°° or p ^ L°° in theorem \3.1l we have 

just to assume that ^r^'VA(p) G L^{L'^). This would imply that p is in fact a weight 
of Muckenhoupt. 

Remark 6 In the sequel the notation of space follows those by Runst, Sickel in 118^ . 

Proof of Theorem I3.lt Our goal is to get a gain of derivative on the density by 
using energy inequalities and by taking advantage of the term of capillarity. We need to 
localize the argument to control the low frequencies. Let ip G C^{W^), we have then by 
multiplying the momentum equation and applying the operator div ( where we use the 
classical summation index): 

dtdiv{ippu) + dij{ippuiUj) - dij{2ipp{p)Duij) - A{ipX{p)divu) + A{LpP{p)) 
= A {AiA{^B{p)) - Al^\VA{p)\'') - Bo,dlj{ipd,A{p)d^A{p)) + 

with: 

d 

-^v = ol^P"^ ' ^"^^ + {di,jf)puiUj + 2di(fdj{puiUj) - {dij(p)p{p)Duij - Al^A{B{p)A(p) 
- 2diipdj{p{p)Duij) - \{p)Aipdivu - 2Vip ■ V(A(p)divu) + Aipap'^ + 2aVip ■ V{p'^) 
- A^{AiAB{p) - Al\VA{p)\^) - • V{AlAB{p) - Al\VA{p)\^) 
+ B^{dlj^)d,A{p)djA{p) + 2Kdi^dj{diA{p)djA{p)) - 2AiA{V^ ■ VB{p)). 

We can apply to the momentum equation the operator A(A)~^ in order to make appear 
a term in AB{p) coming from the capillarity. Then we obtain: 

AlA{ipB{p)) + AlA-\ip\VA{p)\^) + B^A-^RiRj{^diA{p)djA{p)) 
d 

= -A-^—div{ippu) + A-^RiRj{^ppuiUj) - A-\^pX{p)dwu) + A-'^{^pP{p)) (3.10) 

- A-^RiRj{2p{p)Duij) + A-\A)-^R^, 
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where Ri denotes the classical Riesz operator. We multiply now the previous equality by 
A^^'^(ipB{p)) and we integrate on space and in time: 

AiT [ \A'+HipB{p))fdxdt + Al r [ {^\VA{p)\^)A^i^B{p))dxdt 

Jo JRJV Jo JR^ 



+ B^[ [ y R,Rj{^d^A{p)d.jA{p))K'{^B{p))dxdt = 
Jo Jr^Xj 

[ K-^diw{ippu)K^^'{ipB{p)){T)dx - [ A-^div{^poUo)A^+'{^B{po))d: 

Jr!^ JrN 



(3.11) 

J_ f O 

/ iK-^d\Y(ippu)A^+' — i^B{p)) - tfX(p)dwuA%^B(p)))dxdt 
lo Jr^ dt 

+ / / {Ri.j{2ipp{p)Duij)A'{ipB{p)) - RiRj{ippuiUj)A'{ipB{p)))dxdt 
Jo Jr^ 

+ / / ifP{p)A\^B{p))dxdt+ [ [ {A)-^R^A'{ipB{p))dxdt. 
Jo Jr^ Jo Jr'^ 

Now we want to control the term Jq J^n \^^~^^{fB{p))\'^. Before coming in the heart of 
the proof we want to rewrite the inequality ()3.1ip in particular the term: 

A--'div{^pu)A^-^'-^{^B{p)). 
Jr^ '^'^ 

In this goal we recall the renormalized equation for ipB{p): 

^^{vB{p)) + diY{^B{p)u) = -^B{p)diYu + r^, (3.12) 
with r^p = —Vif ■ B{p)u. So by using the renormalized equation (I3.12p we have: 

A-^diy{ippu)A^+' ■^{ipB{p))dxdt = 

A~'^div{ippu)A''{ipB{p)divu) - A~'^d\v{ippu)A'' d\v{(pB{p)u)dxdt (3.13) 

iR^ 

+ / / A^^d\v {Lppu)r^pdxdt. 
Jo Jr'^ 

By combining ()3.1ip and (I3.13P we get: 

Ai r [ \A'+HipB{p))\'dxdt + Al r [ {ip\VA{p)\')A^{ipB{p))dxdt 
Jo Jr'^ Jo Jr^ 



+ Ba [ [ y RiRj{^diA{p)djA{p))A'{ipB{p))dxdt 
Jo Jr«Xj 
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[ A~^div{^pu)A^+'{ipB{p)){T)dx - [ A-^dW{ippouo)A^+'{ipB{po))dx 
Jrn Jrn 

[A""^ div {if pu) A" (if B{p)divu) + A~^div{(ppu)A''div{ipB{p)u))dxdt 

[Rij{2ipp{p)Duij)A'{ipB{p)) - RiRj{ippuiUj)A'{ipB{p)))dxdt (3.14) 

+ / / (pXip)divuA'{ipB{p))dxdt+ [ [ ipP{p)A\ipB{p))dxdt 
Jo Jrn Jo Jr!^ 

[ {A)-^R^A^+i{^B{p))dxdt+ [ [ A-^div{ippu)A^+'r^dxdt. 
Jrn Jo Jr'^ 



In order to control /q J^n \A^^ 2 (^ipB{p))\'^ , it suffices to bound all the other terms of 

2 



(j3.14p . It will allow us to get a control on A^^2 [ipB{p)) and so a gain of | derivative on 



the gradient of density VA{p). 
We start with the most complicated term which requires a control of ^ in L°° if a > — 2 
and a control of p G if a < — 2. In the sequel, we will treat only the case a > —2, the 
proof of the other case follows exactly the same lines. So the fact that ^ is in L°° will 
imply in particular that Vu £ L'^{L'^) and u € L°°(L^), this fact will be permanently use 
in the sequel. 

1) Jo Imi^\^Mp)\'Wi^B{p)): 

By induction we have V{(^B{p)) G L|,(i^2) and by using Sobolev embedding we get 
V{ipB{p)) G LP'{U') with ^ = I — f (we remark that the case s = 2 is critical for Sobolev 
embedding). Now we have: 

^ ^'^^ {2 + a)A{p) (2 + a) ^ ^ 

and we recall that by hypothesis ^ G L°°, so we have ipSJ A{p) G L'^{U') because 
%+a)A[p) ^^ ^ L°°(L'^) for all 1 < r < +00 as V A{p) G L°°{L'^). 

We now consider A* ((/9i?(p)). We have by induction A* ((^i?(p)) G Lp'{H^~i) and A'^{ipB{p)) 
is in L?'{L?') because ipB{p) G L?'{L?') (here the fact to localize by is crucial) which 
enables us to control the low frequencies of A** ((/jS (/?)). We have then A'^{(pB{p)) G 
L2(ifi-f ) ^ L'^{Li) with 1 = 1- Finally by Holder inequality we get ip\VA{p)\'^A%(pB{p)) G 

L^(Li(K^)) because i + J + J = i + i- i + f = land we get more precisely: 

Jip\VAip)\'')A'i^Bip))dxdt < ||VA(p)|U5o(i2)||A^(^S(p))||^2 (^,) 

x||^VA(p)|||^2(^,), (3.15) 

]-u^^l^)\\A'+H^B{p))\\l2^^l,^\\VA{p)U^^^^^^^ 

We proceed similarly for the term: 
rT r 

y^RiRj{ipdiA{p)djA{p))A'{ipB{p))dxdt, 

JR^Xj 



<llil 
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indeed we have in following the same lines ipdiA{p)djA{p) G L^{L'i) with J = 1 - | and 
we use the fact that the Riesz operator is continuous from in for 1 < p < +00. 
We next study the term J^^ A~^div{ippu)A^~^'^ {(pB{p)){t)dxdt. 

2) f^j^ A-^dw{^pu)A^+'{ipB{p))dx: 

We rewrite the term f^^ A~^dw{ippu)A^~^^ {ipB{p)) on the following form: 



A-^dW{ippu)A^+'iipB{p))dx = / A-^div{^pu)A-^+'{ipB{p))dx 

Jrn 

= V / Riiippui)A-'+%ipB{p))dx. 

As i G L^iL"^) then we have u G Lf{L'^). We recall that V A{p) G L°°{H^) then 

ipp G L'^{U') for all 1 < p < +00. We deduce that '^pu belongs to L°°{L'^~'^ n L^) for 
P > 0. So we have Ri{ippui) G L^{L^) for all 1 < r < 2 by continuity of the operator Ri 
from to when 1 < p < +00. 

Case 1 < s < 2: 

Next we have: 

V{^B{p)) = ^B'{p)Vp + B{p)V^ 

then we get V{ipB{p)) G L°°(L2-/5)^ by using the fact that VA{p) G L'^{L'^) and Sobolev 
embedding with Holder inequalities. We have then that LpB{p) belongs to L°°(VF2^_^).So 
A'-^{ipB{p)) belongs to L°°(i/2l|). By Sobolev embedding A*-i(^B(p)) is in L°°{U') 
with i = 2^/3 ~ = 2^ + f ~ ^ with /3 small enough to avoid critical embedding. Finally 
we get Ri{ippu)A-^+'{ipB{p)) G L\,{L'^ {MJ^ ) . Indeed we have | + 2^ = 2^ + f- l<l 
by choosing (3 small enough and ^ + = — 1 + 2^ + f > 1 by choosing f3 small big 
if necessary, we conclude by interpolation. Here the fact that is in a compact support 
is crucial. We have finally: 

A-^d\Y[ippu)A^+'{ipB{p))dx\ <Mo, 
with Mo depending only on the initial data. 
Case < s < 1: 

In this case we conclude by interpolation with the previous case. We now want to study 
the other terms coming from the renormalized equation (|3.12p . 

3) Jq /jgjv A^'^div^Lfp u)A'^ {LpB{p)divu)dxdt, Jj^at A"^div((/?pn)A'*(div((/?i3(p) u))dxdt: 
We start with: 

/ / A-Miv((^/3n)A"-i(div((^S(/3)'u)) = / / d\Y{^pu)A'~^{d\w{^B{p)u)). 
Jo Jm^ Jo Jr'^ 
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Case 1 < s < 2: 
We have: 

div{{ppu) = u ■ V{ipp) + (fpdivu. 

By Holder inequalities and Sobolev embedding we get that di'v{ippu) belongs to L|,(L^~^) 
for all /3 g]0, 1]. Next we rewrite div{ipB{p) u) on the form: 

dYv{ifB(p) u) = u - V{ipB{p)) + ipB{p)diyu. 

As previously div{ipB{p) u) is in L'^{L'^~/^) for all /? g]0, 1]. Now by Sobolev embedding 

we have A*~^div((/3S(p) u) G L\{LP) with ^ = — with (3 small enough to avoid 
critical Sobolev embedding. We conclude that div((/?pn)A*^^(div((^i?(p) u)) is in L^(L^) 
because ^ ~ 2^ ~ ~ 2^ "-'- + §<-'- with (3 small enough if necessary and 

1 + I > 1, so we obtain the result by interpolation. Finally we have: 

i-T r 

A-^dw{ippu)A%dw{(pB{p) u))dxdt\ < Mq, 

Jr^ 

with Mq depending only on the initial data. 
Case < s < 1: 

We have the result by interpolation with the previous case. Next we proceed similarly 
for: 

/ A ^div{(ppu)A'' ^{(pB{p)divu)dxdt. 
/o JR'^ 

4) Last terms 

We now want to concentrate us on the following term: 

/ y RiRj{ippuiUj)A^{ipB{p))dxdt. 
Jm^ 

We know that u G L°°(L2) (as ^ G L°°) and Du e L'^{L^) then u e L^H^) and by 
Holder inequalities and Sobolev embedding we can show that ippuiUj G L'^{L'^~'I^) with 

/? > and so RiRj{(ppuiUj) £ L|.(L^~^). 

We have seen that A%ipB{p)) G L^H^-i) then we have as 1 — | > for /? small enough: 
||A^(^5(p))||^^^^^^<Mo + ||^5(p)||^,^^,_,^ 

with < a < 1. We have then: 

I /" / V RiRj{<ppuiUj)A'{ipB{p))dxdt\ < Mo + \\A^+H^B{p))\\1,,^,. 
Jo JR^ 77 ^ ' 

with < /3 < 1 and Mq depending only on the initial data. 

We are interested in the term: (pX{p)divu A'^{(fiB{p))dxdt We have then divu G 
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LF'{LP') and we have shown that A^{(pB{p)) G L'^{H^~^) so we conclude m the same way 
than the previous term. 

We proceed similarly for the term: Jj^jv Ri,j{'^'^l^{p)Duij) A'^{(pB(p))dxdt. We finally 
conclude with the term: ^ 

/ / ipP{p)k'{^B{p))dxdt. 
Jo Jr^ 

Similarly we have A^{ipB{p)) £ L'^{LP') and (pP{p) G L'^{LP') because ^A{p) is in 
L°°{H^), and we conclude by Sobolev embedding. To finish we have to control the 
term depending on i?^ and r<^ that we leave to the reader. Indeed these terms are easy 
to treat because they are more regular than the previous terms. 
We finally get by using all the previous inequalities: 

||Ai+i(v?i?(p))||i.(^.) < Co||V^(p)||ioo(i2)||Ai+f 

+ Ci||Ai+t((^i?(p))||2^(^,) + Mo 

with < /3 < 1 and Co, Ci, Mq depends only of the norm of initial data. By energy 
inequalities we have ||V^(p)||ioo(2,2) < e < 1, we can then conclude that: 

\\^B{p)\\^,^^,+^^<Mo 

with Mq depending only on the initial data and ip. □ 
We want now improve this result by extracting a specific structure of the capillarity 
tensor. Indeed by choosing n{p) = Hip~^ with k > 0, we are going show that we get the 
same estimate but without any conditions on the vacuum. In fact, the power —2 seems to 
be critical, and it explains very easily in our proof by the simple fact that A{p) = B{p). 
So we can apply a boobstrap argument without asking a control of p or ^ in L°°. 
However a new difficulty appears on the control of the pressure term as in our case we have 
only a control on V Inp G L°°{L?'). So we can hope directly a control of tpp G L°°{LP) for 
all 1 < p < +00. We will see in fact that we can in the same time as a gain of derivative 
on the density obtaining a gain of integrability on the density. 

Theorem 3.2 Let N = 2, k > and {p,u) be a smooth approximate solution of the 
system ip with k{p) = Kp^^ .Then there exists a constant r/ > depending only on the 
constant coming from the Sobolev embedding such that if: 

I|Vpo||l2(R2) + ||^|lto|||L2(]g2) + ||i7(po)||Li < V 

then it exists a > such that for all tp G C^(M^); 

+ yp"~^^p\\hiL^) + ll'/'PllIt+.((o,T)xRiV) < M With 0<s<e, 

where M depends only on the initial conditions data, on T , on p and on e. e depends 
only of 7 the coefficient of the pressure and is small. 

Proof: The proof follow the same line as the proof of theorem 13.11 except concerning 
the bounds of estimate coming from the capillarity term and the fact that we lost the 
control of p in L^{L^^^) for all 1 < p < +oo. We need then to get a gain of integrability 
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on the density to treat the term coming of the pressure. We apply to equation ()3.9p the 
operator A and we multiply par p° with a > 0: 

[AiA{ipB{p)) + Al{ip\VA{p)\') + B^R,Rj{ipd,A{p)d,A{p))] p" 

= ^^{A)-^d\w{ippu)p'' + R,Rj{^puiUj)p'' - v?A(/5)(divn)p" + (3.16) 

- RiRj{2p{p)Dui^j)p^ + (A)-ii?^p", 

Next we integrate on (0,T) x the equation 13 . 1 6 1 we get then: 
I ^P^^"{x, t)dxdt + A^f p'^~^\Vp\^dxdt = 



'(0,T)xRJV J{0,T)xM^ 

/ 

'(0,T)xRJV 



Ba / ipdiln{p)djln{p)Rijp°'dxdt - / {A)-^dw{ippu)p°'{T,x)dx 

+ / {A)-^div{ippouo)podx- / {A)-^dw {if pu)dtp'^dxdt (3.17) 



Rij{ippuiUj)p°^dxdt — j (pX{p){divu)p°^dxdt 
- / RiR,j{2p{p)Duij)p°'dxdt + / {A)-^ R^p'^dxdt. 

In the sequel we will note: 

{A)-'^dw{ippu)p''{T,x)dx+ / {A)-^div{ippouo)p^dx 



(A) ^dYv{ippu)dtp"dxdt — l Rij{LppuiUj)p"dxdt 

(0,T)xRJV J(0,T)xRJV 



ipX{p){divu)p'^dxdt — / RiRj{2p{p)Duij)p'^dxdt 

(0,T)xR^ y(0,T)xR^ 



+ / {Ay^R^p'^dxdt. 



'(0,T)xI 

Our goal is now to have p £ Lj^" ({0,T) x M^), for this we have to control all the terms 
on the right handside. This processus follow the same lines that this one used in the case 
of Navier-Stokes compressible to get a gain of integrability on the pressure. The new 
difficulty consists only in the following term: 



(fdi ln{p)dj \n[p)Rijp^dxdt, 

(0,T)xR^ 
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To control this term we need of a gain of derivative on Vln(/9) that is why we have then 
in suming equation (|3.17|) and (j3.1ip : 



All I \^^^HvB{p))fdxdt+ [ {ipp^+"{t,x) + Aap'''^\Vp\^)dxdt = 

Jo JR^ J(0,T)xR'v 

-All [ {^\VA{p)\'^)A'{ipB{p))dxdt- Be, [ ipdiln{p)djln{p)Rijp"dxdt 

Jo JR^ J(0,T)xIR^ 

-Bo, [ [ y RiRj{ipdiA{p)djA{p))A'{^Bip))dxdt + F 

Jo iR^ ~f 



A~'^div{^poUo)A^-^'iipB{po))dx + / A~^div{ippu)A^+'{^B{p)){T)dx 

d 



T 



{A~^dw{ippu)A^+' — {^pB{p)) - ipX{p)divuA%ipB{p)))dxdt 

Jr^ 

+ / / {Rij{2ipp{p)Duij)A'{ipB{p)) - RiRj{^puiUj)A'{ipB{p)))dxdt 
Jo Jr^ 

+ [ [ ^P{p)A\ipB{p))dxdt+ [ [ {A)-^R^A'{^B{p))dxdt. 
Jo Jr^ Jo Jr^ 



We have now to control all the term on the right handside. We begin with the capillarity 
term: 

T 

^\V{lnp)\'^A%ip\n{p))dxdt, 

We have then by induction and Sobolev embedding V{iplnp) G L'^{U') and A''((/?ln(p)) 
is in (L^i ) where | = | — | and ^ = | so by Holder inequalities we have a control of 

((5|V(lnp)pA''((^ln(/9)) in L\,{L^) because ^ + f + ^- f = l- Next we treat the following 
term: 

ipdi \n.{p)dj ln{p)Rij{p")dxdt, 

(0,T)xR'v 

We have then by induction Rij{ipp°') G L^{{0,T) x M^) and (pdi{lnp) G L'^{LP) 
with p = |, by Holder inequalities we have ipdiln{p)djln{p)Rijp°' G L^(L^) because 
^ + l + ^<landi + ^<l. 

Similarly we have ipP{p) G L^((0,r) x M^) and Ann{ipp) G L°°{L-s) by Sobolev 
embedding, we conclude by Holder inequalities as + f < 1- To finish we study the 
term asking control L°° in time as A~^di\{ippu)A^'^^^ {Lphi{p)){T)dx coming from 
F, we have: A-^d\Y{ippu) G L°°{LP) and A^+''{Lp\n{p)) G L~(L'?) with ^ = ^ + | 
and i = i - ^ = f - i, we have then A-^d:iY{ippu)A-^+'{^\n{p)){T) G L°°{L^) if 

^ + I < 1. The others terms are left to the reader. We conclude by writing the final 
estimate where we set: 



A{T)= [ ipp^+''{x,t)dxdt + A^ [ p'^-^lVpl^dxdt, 

B{T)= r f \A^+^{^B{p))\^dxdt. 
Jo Jr^ 
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We have then the following final estimate: 

AiT) + B{T) < eB{T) + eA{T)^ B{T)^ + CA{Tf^ + C' B{Tf^ + M, 

with < /3i,/32 < 1- By boosbstrap we can conclude. □ 

Remark 7 The major difficulty in the previous proof is to treat the terms coming from 
the capillarity, it implies to impose a smallness condition on the initial data. An other 
idea would be to use refined Sobolev embedding to avoid smallness condition. So we have: 

I ^\V\n{p)\''Knn{p)dxdt < ||v9Vln(p)||i2(iPi)||AMn(p)||i2(i,i)||Vln(p)||ioo(i2), 

J(0,T)x9.'^ 

with ^ = ^ ~ f one? ^ = |- Next by using propositions \2. 51 we get: 

\\ipVHp)\\L.<C\\^\np\\^s llv^lnpll J^. 

-'-'00,00 

\W^"Hp)\\lv < cy\np\\~\ yinpwi^^s 

-'--'00 ,00 

We have then: 

||v?Vln(/))||iPi ||v?AMn(p)||L9i <C'^\\Lphip\\ s ||(^ In j 

^00 ,00 

We could conclude our argument if ip In p € L'^{B^^) with e arbitraly small. In this 
case we would be able to get a gain of derivative on the density without any condition of 
smallnesse on the initial data. We can notice that in f^, p is in L'^{B^ i) so it's wide 
enough. 

In the following result, we want show that we can obtain similar result when we approx- 
imate the capillarity coefficient k{p) by a constant k. 

Corollary 1 Let N = 2 and a, M, k € M. {p, u) is a smooth approximate solution of the 
system (I j. with the following capillarity coefficient: 

where 61, 62 are regular function such that n is a regular function. Then there exists a 
constant rj > depending only on the constant coming from the Sobolev embedding such 
that if: 

I|Vpo||l2(R2) + ||^|no|||L2(]g2) + ||i7(po)||Li < V 
then we get for all ip G C^(M^); 

11(^5(^)11^2^^^1+1^ < M with0<s<2, 
where M depends only on the initial conditions data, on T , on ip and on s. 
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Proof: The proof follows the same lines than proof of theorem 13.11 We have then in this 
case: 

A{p) = logpl{p<a} + /3l(p)l{Q<p<2a} + I^P^{2a<p}, 
B{p) = logpl|p<„} +/32(p)l{a<p<2a} +/'^l{2a<p}- 

where [ji{x) = V9i{y)l{a<p<2(x}dy and (32{x) = y6'i(?/)l{a<p<2a}(iy.Here to simplify 
we assume that the viscosity coefficient are constant. At the difference with the previous 
proof, this case is more simple because we have ^{pl{2a<p}) ^ L°°{L?'), so we have easily 
p G L°°{Lf^^) for all 1 < p < +00. It implies that we don't need of condition of smallness 
on s as in the previous case. 

We apply now to the equation (j3.9p the operator (A)^^div, next we multiply by A{B{p)) 
and we integrate on space and in time: 

, „A^+*i?(p)P + \VAip)\\Pf\{V{^B{p)) - Bip)VipW{^Bip))dxdt 

Jm B Kp) 

+ [ [ y"RijiKip)^dipdjp)A'{^B{p))dxdt= [ A-^div{ippu)A^+'{^B{p)){T)dx 
Jo Jr^ ~1 ' J^'^ 

''iJ 

A-^dW{^PoUo)A^+'{^B{po))dx - [ [ A-^diw{ippu)A^+' ^{^B{p))dxdt 

Jo Jr'^ 

+ (2/i + A) / / ipd\YuA\LpB{p))dxdt- / / RiRj{(ppuiUj)A%ipB{p))dxdt 

Jo Jr'^ Jo Jr'^ tI 

+ [ [ tfP{p)A'{ipB{p))dxdt+ [ [ {A)-^R^A%(pB{p))dxdt. 
Jo Jr'^ Jo Jr^ 

Now we proceed as in the proof of theorem 13.11 and we have to control the terms: 

, {\VA{p)\{Pf\{V{ipB{p))-Bip)Vip)A%ipBip))dxdt, 
Jr^ B [p) 

Tre crucial point is that is in L^{L°°). So as previously we have: 

i\VAip)\^(V{ipBip))A%^B{p))dxdt 
B (p) 



JR'^ 



<C\\^^\\l^^l^^\\VA{p)\\l^^^l^)\\V{^B{p)^^^^^^^ 

<C||^llL-(L-)l|V^(p)||Lo.(^.)||Al+f(v.5(p)||i.(^.), 



with I = 1 — I and | = |- We get finally: 

f [ \^VA{p)\'A^{^C{p)) < 4l|Ai+^(^i5(p))||2^ 
Jo Jr^ 

For the other term we proceed similarly as theorem 13. 1[ □ 
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3.1.1 What happens when we choose a coefficient of capillarity k large 

Assume that K,{p) = Kp" with a > and k is a very large constant coefficient. We have 
then by energy inequality: 

||V^(p)boo(^2) < C{\\VA{p^)\\l2 + -(||^no||i2 + IIpIiI,)). 

It means that we can hope gain of derivative on the density if we assume that k is enough 
big and only V A{pq) small in L?'. We obtain then the following theorem: 

Theorem 3.3 Let N = 2 and a, M € M. (p, u) is a smooth approximate solution of 
the system with initial data in the energy space and with the following capillarity 

coefficient: 

where 6i , 02 are regular function such that k is a regular function. Then there exists a 
constant r] > and k enough big depending only on the constant coming from the Sobolev 
embedding and on the initial data such that if: 

||VA(/)o)||l2(r2) <ry 

then we get for all ip G C^(M^); 

\\ipB{p)\\^,^^^,+ s^ < M with < s < 2, 

where M depends only on the initial conditions data, on T , on ip and on s. 
Proof: In fact by following the proof of theorem 13. 11 we just have to check that: 

kT f ^\VA{p)\''K\ipB{p))dxdt < K\\Bip)f s 
And this is the case because with our hypothesis: 

\\VA{p)\\l^^l^) < Ci\\VA{po)\\L2 + -i\\VP^uo\\L2 + \\p\\%)) < e, 
with £ enough small when k is enough large and ||V74(/3o)||l2 enough small. □ 

3.2 Control of the density without any conditions of smallness on the 
initial data 

In this section, we want avoid to impose some conditions of smallness on the initial data. 
We would get a theorem of global weak solution with large initial data. For making, we 
have to localize our previous arguments with test function E Cq° with small diameter 
support. Indeed we would get a control of \\lB(x,r)^ ■^{p)\\l^(l'^) aiid prove that this 
norme is small uniformly for all x G when r is small. So for a chosen compact K, we 
have to split it in a finite union of small ball and apply our previous argument on each of 
these balls. Before entering in the heart of subject, we would start now with localizing 
the classical energy inequality on small support. 
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3.2.1 Localizing energy inequality and phenomena of concentration 

We are interesting in getting new energy inequalities which take in account the support of 
the initial data. In particular we want investigate what happend when we choose initial 
data localize in small ball. So for making we let ip G C^(R^) with (/? = 1 on B{xo,R) 
and supp(/7 C B{xq,2R). Multiplying the equation of momentum conservation in the 
system (jl.ip by ipu and integrating by parts over (0,t) X R^, we obtain the following 
estimate: 

I I ^{x){p\u\\t,x) + ^\Vp\'' + 2Il{p) - 21\{-p)){t,x)dx 

+ I I (/^|Vn|^ + ,^|divn|^)(^d3; < / ^(f{x){poul + k{po)\V po\'^)dx 
Jo JR^ JR^ ^ 

+ / (/?(x)(n(po) -n(p))(x)d2;+ \ [ \{u-Vu).Vip + P{p)u-V^]dxdt (3-18) 
Jrn Jo Jrn 

+ K{p)Vp-Vipdw{pu)dxdt+ / / {k{p) + ^pK'{p))\Vp\'^u-Vipdxdt 

Jr Jr'^ Jr Jr^ 2 

+ / / pK{p)dj pdjUidiip dxdt + / / pK{p)djpuidijip,dxdt. 
Jr Jr'^ Jr Jr^ 

Now we have to control the term on the right handside in the goal to get some energy 
estimate localized in space. We recall that in the proof of theorem l3.H we need a condition 
of smallness on ipVA{p) G L°°{L?'). So by a condition of small support on ip we can hope 
get this condition of smallness on ip\I A{p) G L°°{L?') for a finite time. In fact we search to 
prove that there is not concentration effect in some point of the space for a small intervall 
[0,T]. More precisely it means that for any sequel of approximate solution {pn,Un)neN, 
we have the following property: 

yK a compact, 3e > Vt G [0,T] \\V {A{pn){t, ■))1k{-)\\l2 < e. 

So we need of this type of property to can use a boobstrap in the proof of theorem 
13.11 without asking condition of smallness on the initial data. Moreover this property is 
very natural. Indeed we want avoid some concentration effect. In particular if (pn,u„) 
is a sequel of approximate solution for the system (II. ip . we want ask that |VA(/9„)p 
doesn'tconverge to some Dirac measures. 

Proposition 3.7 Let q > 0, /3 > Q,r] > 0, k > 0, a regular function such that: 



Let (p, u) a regular approximate solution of system with large initial data in the 

energy space. Assuming that u G L}p^^{L°°) n L^+'^(L^) then it exist a time Tq > such 
that for all ip G C^(M^) we have: 

where M depends only of the initial data, K and ||u|Ll+/3/roo^• 
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Proof: The proof consists only to get a control of ||V(A(/9)) li^||2,oo(-2,2) for K a compact 
of small enough measure. Indeed by following the same lines as the proof of theorem 13 .11 
the main difficulty is to control the term for ip G C^(M^): 

/ ip\VA{p)\'^A'B{p)dxdt. 
Jr^ 

In our case, we have to take in consideration that the integration is only on the support 
of (p. So let K = suppip a compact of M^, we consider a subordinate partition of unity 
of the indicatrice Ik with {(pk)i<k<p of class C°° such that: 

• supp(/5fe c B{xk, A) C K, Ylk 95fc = 1 on K and < v9a: < 1, 

• ||V"(^fc||Lcx, < C„|A|-° for \a\ < 2, 

• llVv^fclliP < Ca ] 1 for 1 < p < oo. 

We have now to estimate the following term: 

rT r 

^k\VA{p)\^K'{ipkB{p))dxdt = 

JR^ 

f [ VA{p) ■ {^,^VB{p))A%^kBip))dxdt + R^, 
Jo Jr^ B [p) 



(3.19) 



to control ipkB(p) in Lp'i^H^^^) and where here Rip^^ is a terme easy to control. We want 
use a argument of boobstrap as in the proof of theorem 13.11 We have now just to prove 
that for e enough small, it exists a Tq such that we have: 

II lsupp<^feVA(p) 11^00(^2) < e. 
Let V'fc S such that V'fc = 1 on suppi^fc and suppV'fc C 2supp99fc We set then: 
A{T,^k) = \ I Mx){p\u\\t,x) + \VA{p)\^ + 2Uip)-2U{p)){t,x)d 
By the equation (|3.18|) . we have: 



X. 



A{T,^Pk)+ I I {p\Vu\^ + ^\divu\'^)^kdx < I l^lJk{x){poul + \VA{po)\^)dx 
10 Jr^ ■'^'^ ^ 



+ [ MxMipo) -'n{p)){x)dx+ [ [ [{u-Vu).Vi:k + P{p)u-VyJk\dxdt 

JRN Jo JR'^ 

I I K,{p)V p ■Vil)kd\Y{pu) dx dt + I I [k{p) + -pK {p)^\V pI'^u ■ Vipkdxdt 
Jr Jr'^ Jr Jr^ 2 



+ 

+ 11 pK{p)djpdjUidiTl>kdxdt + I I pK{p)djpuidi,jilJk,dxdt. 
Jr Jr'^ Jr Jr^ 

We have to control all the term of right to have estimates on A(T,'ijjk). Let e > 0, easily 
for A enough big, we have: 

/ l:Mx){poul + \VAipo)\^)dx+ [ Mx)(n-{po)-U{p)){x)dx<^. (3.20) 
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We have now to control the other terms, we begin by the most compUcated and we have: 

/ / Hc{p)Vp ■Vipk'iiv{pu)dxdt < \\^A{p)\\l^,^2)\\Vi;k\\L°°\\u\\L^+i3(L'^)T^ 
JrJrn ^ ' 

+ ||V^(/5)||ioo(i2)||VV'fc||L-||Vn||i2(i2)T3. 
We can treat similarly the other terms and we get finally by using inequality (j3.20p : 

^2-/3 \\yu\\L'2{L'^)T'2 

+ ll"llL2(/fl)||M||L2+/3(L2)||VV'fc||L°°r2(2+/3) + 1 1 P(p) 1 1 ^oo (^1 ) 1 1 U 1 1 ^1+0 (^oo) 1 1 VV'fc 1 1 ^ l+Z' , 

By the conditions on and the previous inequality we get: 

A{T, V'fc) < I + C|A| (||u||ii+,(ioo)T^ + + \\u\\l2+0^l^^T^) 

+ \Ml^+0(l-o)T^ + |A|r5). 

For a small enough time Tq depending of the initial data, ||^i||Li+ff(Loo'), ||ii||L2+/3(j;^2) and 
A we have: 

A{T,iPk)<e. (3.21) 

We can now come back to the crucial point of the proof of theorem I3.H and from the 
equality ()3.19p and inequality (I3.2ip we have: 

lo Vk\^A{p)\^A%ipkB{p))dxdt < CUk^A{p)\\Lo.^L^)ykB{p)\\l,^^,+ s^ + Mo, 

< CeUk^A{p)h^^^2)ykB{p)\\l,^^,^s^ + Mo. 

By choosing e enough small, we conclude that ipkB{p) E L^(i?^^2). Next we get by 
Sobolev embedding ipkVB{p) G L'^^{LP) with | = ^ - |- So we have Ik^B{p) = 
X]fc fk^-^ip) which is in L?'{U'), by coming back at the proof of theorem 13.11 we finally 
get a control of ^B{p) in L'^{H^^^ and this conclude the proof. □ 

Remark 8 In a similar way, we could show that if we control the high frequencies of 
VA{p) in L°^{L?') then we can get a gain of derivative with large initial data in the energy 
space. More precisely we have: 

diA{p) = diA{p)l{i^i<M} + 5i^(p)l{|5|>M} = /i + /2- 

We know that fi is regular and in particular fi G L°°{LP) with p G [2,+oo[, so by 
using the same argument as in proff Wl\ we are able to treat this term. The main 
difficulty is to control 9jyl(p)l{|g|>j\,f}, indeed we would use a argument of smallness on 

||5i^(/3)l{|g|>M}||L°°(L2) for M enough big. 

In particular for an approximate sequence of solutions, is it possible to find for a e > 
M enough big and depending only on the initial data such that: 

||<9j^(/3„)l{|^|>7V/}||L°°(L2) < ^• 

The difficulty is to show that the mass does not run away to the inifinity. 
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3.3 Existence of global weak solutions for N = 2 

We may now turn to our compactness result. First, we assume that a sequence {pn, Un)n<m 
of approximate weak solutions has been constructed by a mollifying process, which have 
suitable regularity to justify the formal estimates like the classical energy estimates and 
these coming from your previous theorems. 

Remark 9 In fact, we easily construct this sequence of solution {pn,Un)nem by using the 

JV — — 1 

result of jy^, it suffices to choose initial data {p^,u^) G x with a condition of 

smallnees. We know indeed that there exists tj > such that if: 

WpIW n + hnW N < r], 

then the solution [p^,u^) are global and strong. 

Moreover in the sequel the viscosity coefficients check the properties of definition 13.51 it 
means that they are non degenerate. Moreover in the sequel we will work in dimension 
2, and to simplify we choose a isentropic pressure P{p) = p^ with 7 > 1. 
Moreover this sequence {pn,Un)neN has initial data i{po)n, iuo)n)) close to the energy 
space. By using the above energy inequalities, we assume that j'y{{po)n), I ^(^((/'o)™) P 
and (/5o)n|(^o)nP are bounded in L^(M^) so that (po)n is bounded in L^i^^) ■ 

Then it follows from the energy inequality that: 

1- hipn), |VA(p„)|2, p„|n„p are bounded uniformly in L°°{0,T, L'^{R^)), 

2. Dun is bounded uniformly in L'^{R^ x (0,r)), 

3. Un is bounded uniformly in L'^{0,T, (Br)) for all R,T £ (0, +00). 

Extracting subsequences if necessary, we may assume that pn, Un converge weakly re- 
spectively in L°°(L^(R^)), L'^{0,T;H^{Br)) to p, u for all R,T £ (0,+oo). In a similar 
way to the case of Navier-Stokes compressible, we can extract subsequences such that 
y/p^Un, PnUn, PnUn ® Un Converge weakly to ^/pu, pu and pu ® u. 

In fact the main difficulty is to verify that the quadratic gradient term of the density 
VA{pn) ® VA{pn) converge to V A{p) ® VA{p). 

Finally we add uniform hypothesis on B{pn) coming from the gain of regularity that we 
have obtained in the previous part, so we assume that: 

V(/7GC^(M^), is bounded uniformly in with s > 0. 

Moreover when Hi{p) = with k > we have: 

• ippn is bounded uniformly in Li'"^°((0,T) x M^). 

• ippn~'^^ Pn is bounded uniformly in L^(L^(R^)). 

We can now show the two following theorem. The next theorem concerns the existence 
of global weak solutions in the case where k{p) = Kp~'^. 
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Theorem 3.4 Here we assume k{p) = with k > 0. There exists r] > such that if: 

||V(ln;.^)||^. + Wy^Klh^ + \\jM)\\L^ < r/ 

then, up to an extraction a subsequence {pn,Un) converges strongly to a weak solution 
{p,u) (see definition \2.4\ ) of the system il.l]) . Moreover we have Vlog(/9„) (Ei Vlog(p„.) 
converges strongly in Lj^^{M. x M^). In addition p check for all tp G C^; 

where M depends only on the initial conditions data, on T, on ip and on s. 

Remark 10 This theorem is a theorem of existence of global weak solution and not only 
a result of stability. Indeed as explained in the remarka O we are able by using the result 
of to construct global approximate solution. And the result of ^9] is compatible with 
our hypothesis of smallness. It would be the same for the following theorems. 

The next theorem treat of global weak solution for capillarity coefficients approximating 
a constant. 

Theorem 3.5 Let N = 2, a>0, e>0 and the following capillary coefficient: 

where 6i , 62 are regular function ssuch that k is a regular function. There exists i] > 
such that if: 

l|v(A(po"))||^. + \\y^K\h2 + \\jM)\\L^ < V 

then, up to a subsequence {pn,Un) converges strongly to a weak solution {p,u) (see defi- 
nition {2^ ) of the system Moreover we have VA{pn) CS> VA{pn) converges strongly 

in LI^^{M. X R^). In addition p check for all ip G C^; 

\\'^^^P^\%{n^^^) - ^'^^ < s < 2, 

where M depends only on the initial conditions data, on T, on ip and on s. 

To finish, we give a theorem of existence of global weak solution when the capillarity 
coefficient k is big and when only Vj4(/9o) admits a condition of smallness. 

Theorem 3.6 Let A^ = 2, a>0, e>0 and the following capillary coefficient: 

^ ^^2+F''^{P<"} + ^l(/')l{a<p<2o} + '^l{p>2a}- 

where 61, 62 are regular function ssuch that k is a regular function. Let ^/pqUq G L'^, 
j-yi^Po G L^ There exists k enough big and r] > such that if: 

\\viAip^mL^<v 

then, up to a subsequence {pn,Un) converges strongly to a weak solution {p,u) (see defi- 
nition \2^ of the system Moreover we have 'VA{pn) VA{pn) converges strongly 

in LI^^(R X M^). In addition p check for all p G C^: 

\\'^^^pKuh'+^)-^^ mi/i0<s<2, 
where M depends only on the initial conditions data, on T, on ip and on s. 
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Proof of the theorem 13.41 : The main difficulty states in proving the convergence of 
the following nonlinear terms P{pn) and V Inp^ V ln/?„. The other terms follows the 
same lines as the compressible Navier-Stokes problem studied by P-L. Lions, E. Feireisl, 
A. Novoty and many other authors. 

According to theorem 13.11 we have seen that for all 99 e C^(M^) iflnpn G L^H^+I). 
So we have Vv? e C^(R^) 

We can now use some results of compactness to show that V(lnp„) converge strongly in 
L'^{Lf^^) to Vlnp. We recall the following theorem from Aubin-Lions ( see Simon for 
general results [19j). 

Lemma 1 Let X ^ Y ^ Z be Hilhert spaces such that the embedding from X in Y 
is compact. Let {fn)neN a sequence bounded in L'^{0,T;Y), (with 1 < q < +00) and 
bounded in LP{0,T;Z) (with 1 < p < +00), then {fn)n&N is relatively compact 
in Li{0,T;Y). 

We need now to localize our arguments because we want use some result of compactness 

s_ s_ 

for the local Sobolev space Hf^^^ with s > 0. H^^^^ is compactly embedded in Lf^^. Let 
(Xp)pGN be a sequence of Co°(M^) cut-off functions supported in the ball B{0,p + 1) of 
and equal to 1 in a neighborhood of B{0,p). We have then by using mass equation: 

^V(lnp„) + Vdiv(u„) + V('u„- Vlnp„) = 

We can then show that (^(xpV(ln/9„)))^gp^ is uniformly bounded for all p in L'^{H°') for 
a < by using energy inequalities. Moreover (xpV(lnp„))^^j^ is uniformly bounded for 
all p in L^(//2). Applying lemma [D with the family (xpV(ln/9„))^gj^ and X = XpH^, 
Y = XpL'^j Z = XpH"' and using Cantor's diagonal process, we provides that after up to 
a subsequence: 

Vp > XpV(ln/>„) ^„^+oo Xpa in L^(L^). (3.22) 
with for all p G N , XpO- G L'^{Hi). 

Moreover as \yf]h(S Pn = yfPn S L^{LP) with ^ = ^ + by the same argument we 

have ^/p^ converge strongly to a certain b in L^(L^^~^) with e > such that 27 — e = 2. 
Moreover up a subsequence ^/fhi converges a.e to b. We have then for all ip G C^: 

/ (pp'^dxdt — >„^_|_oo / ipb^dxdt. 

And as pn converges weakly to p, we have shown that b^ = p and that pn converges a.e 
to p. 

We can show that for all (p G ip In p„ converges weakly to iplnp by the fact that pn 
converges a.e to p and the proposition 12.61 It means that a = Vlnp. Finally we have 
shown that V(ln/9„) converges strongly to Vlnp in L^(L^^^). We have then obtained 
that V In Pn (g) V In pn converges in distribution sense to V In p (g) V In p. 
The last difficulty is to treat the term P{p"'), we proceed similarly. Let if G C^(K^), 
so as ip\J yfp^ = ^ip^/p^V In pn and as p„, is uniformly bounded in L'p'" {Lj^") we get 

Vy/p^ is uniformly bounded in L^'^~^'^\l^^J with ^ = | + 2('r+a) • ^^^-^ conclude by 



p 2 ^ 2(7+0) • 

compact Sobolev embedding and proposition 12.61 □ 
The proof of theorem 13.51 and 13.61 follows the same line than the previous proof. 
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4 Existence of weak solution in the case = 1 

We are now interested by the case = 1. To start with, we focus on the gain of derivative 
for 'VA{p) with a general capillarity term. 

4.1 Gain of derivative 

We can now write a theorem where we reach a gain of derivative on the density p by 
using the same type of inequalities as in the case N = 2 and without any conditions of 
smallness on the initial data. 

Theorem 4.7 Let {p, u) he a regular solution of the system M.l\) with initial data in the 
energy space and let n{p) = Kp'^ a general capillarity coefficient with a G M. Then we 
have: 

with < s < ^ and Mq depending only of the initial data. 

Remark 11 We observe the two important facts: 

1. We don't need any hypothesis on the size of the initial data. 

2. We don't need to localize because we know that p G Lf'^. 

3. We don't need to assume that ^ G L°° . 
Proof of theorem 14.71 : 

We use here the same estimates as in the previous proof except for the delicate term: 
lo Ir \^xAip)\'^^'^B{p). In the sequel we will show only the case a > —2, the proof of 
the other case is similar. We have then dxA{p) G L°°{L?') and p — p ^ L°°{Lp') so by 

Sobolev embedding A{p) G Next we set VB{p) = A{p) = \A{p)VA{p). 

But we know that A{p) belongs to L°°(L°°), so that VB{p) G L°°(L2). Finally we get for 
< s < 1, A*5(p) G L°°{H^-'). Now for < s < i by Sobolev embedding we obtain: 

A'B{p) G L°°(L°°). 

So we can control the term Jj^ \dxA(p)\'^\A'^ B{p) as follows: 

r [ \dxA{p)\'\A^B{p)\ < \\dxA{p)\\i^^^,y 
We treat the other terms similarly as in the previous proof. □ 
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4.2 Results of compactness 

We can now prove our result of stability of solution in the case = 1 by using the 
previous gain of derivative. Let {pn,Un)nm a sequel of approximate weak solutions of 
system (II. ip . We get now a theorem of global weak solution with large initial data. 

Theorem 4.8 Let (/9q,Mq) initial data of the system 1^1. 1\) in the energy space Let k{p) = 
~2TFl{p<a} + ^{p)^{a<p<2a} + ^l{<p>2o} ^^^^^ ^ > and 6 a regular function. 
Then up to a subsequence, {pn, Un) converges strongly to a weak solution (p, u) on (0, T) x 
M in the sense of the distribution for all T S (0, +cxd) (see definition \2.4^ . Moreover 
dxA{pn) converges strongly in L^(0, T, L^^^(M^)) to dxA{p). 

Proof:The proof follows the same lines as in the proof of theorem 13.41 □ 



5 Appendix 

5.1 Computation of the capillary term 

This section is devoted to rewrite clearly the capillarity tensor K in the the goal to express 
the non linear terms in distribution sense. On other interesting reason is to describe the 
regularizing part of the capillarity, it will allow us to extract smooting effect and so to 
treat the non linear terms in distribution sense. We recall that: 

divK = V{pk{p)Ap + ^{k{p) + pk'{p))\Vp\^) - div {k{p)V p ^ V p) . (5.23) 
and as for all / € C^: 

Af{p) = /{p)Ap + f"{p)\Vpf. 

We get then: 

divK = V{Af{p) - ^{k{p) + p>,'{p))\Vp\^) - div{K{p)Vp ® Vp). 
with f'{x) = xn{x). It gives in particular the estimates (jl.3p . 

5.2 Inequality energy estimates 

We are interested here in derivating bounds estimates on the system (II. ip . We have 
to multiply momentum equation by u and integrate over the time and the space. We 
concentrate us only one term: 

/ / divK ■ udxdt = / / ( pK{p)diAp + pK {p)dipAp + ^pk" {p)dip\V pf 
Jm. Jr'^ Jr Jr^ V 2 

+ pK {p)djpdijp^ Ui dx dt 

Next we have: 

(^/"^"(p)f^iP|VpP + Pi^ {p)djPdijp)ui dxdt = 

- K {p)\Vp\^u -Vpdxdt — - / / pK {p)\Vp\'^diYudxdt, 

2 Jr Jrn 2 J]g J^N 
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and: 



/ / [pK{p)diAp + pK {p)Apdip)uidxdt = 

~ K{p)Apu -Vpdxdt — / / pK{p)ApdiYudx dt. 

Jr jr^ Jr Jr^ 

By mass equation we have: 

pdivu + u ■ Vp = —dtp. (5.26) 
In using ([EMI), <^^^ and (lOHl) we get finally: 

/ / divK-udxdt= f f K{p)Apdtpdx dt + f f k {p)\Vp\'^dtpdxdt, 
Jr jr^ Jr Jr^ 2 J^n 

/ / K{p)\Vp\'^ dxdt. 
Jr Jrn 



ld_ 

''2dt 



The other terms are classical. 
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